Experimental quantum-enhanced estimation of a lossy phase shift 
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A paradigm for quantum-enhanced precision 
metrology is found in optical interferometry 
capable to sense physical quantities as diverse as 
position, time delay or temperature through a 
measurement of a phase shift. When standard 
light sources are employed, the precision of the 
phase determination is limited by the shot noise, 
whose origin can be traced to the random man- 
ner in which individual photons emerge from the 
interferometer. Quantum entanglement provides 
means to exceed this limit [2], |3|, y, |5|,[6| with the 
celebrated example of NOON states 0, 0, I3> E3| 
that saturate the ultimate Heisenberg limit on 
precision [11], but at the s ame time are extremely 
fragile to losses [12;, Il3l Il4l |. In contrast, we 
provide experimental evidence that appropriately 
engineered quantum states |l5j | outperform both 
standard and NOON states in the precision of 
phase estimation when losses are present. This 
shows that the quantum enhancement of metrol- 
ogy is possible even when decoherence is present, 
and that the strategy for realising the enhance- 
ment is quite distinct from protecting quantum 
information encoded in light [16l. Il7j |. 

The Mach-Zehnder interferometer shown in Fig. [TJa) 
serves to illustrate the principle of phase estimation. The 
input light is divided into two beams and one of them 
probes the phase shifting element. Recombining the two 
beams maps the acquired phase information onto the 
output intensity, read out via photodetection. The pho- 
todetection signal exhibits fluctuations which result in 
the shot noise limiting the measurement precision. In 
the semiclassical picture this effect is attributed to the 
uncertainty in the photoelectron number generate d by 
an incident beam even if its intensity is constant [18j |. 
In the fully quantum picture, light is composed of ele- 
mentary quanta — photons — whose behaviour can be pre- 
dicted only statistically, hence leading to fluctuations at 
the interferometer output Q. The shot noise scales as 
1 / \/N, where N is the average number of photons used 
for the measurement. For a lossy phase shift, it is nec- 
essary to readjust the splitting ratio for the input light 
which yields a generalized bound on precision [1 91 ] in the 
form of the standard interferometric limit (SIL). 

More generally, the two beams — one sent through the 
sample and the second one serving as the reference — can 
be prepared in a certain probe state \ip). According to 




FIG. 1: (a) The generic Mach-Zennder interferometer fed with 
a classical light beam. The phase shift ip inside the interfer- 
ometer modulates the output intensities. Possible losses ac- 
companying the phase shift are represented by a beam splitter 
with power transmission n; (b) Linear network to demonstrate 
optimal phase estimation. The beam splitters with power 
transmission i?i and $2 prepare the optimal probe state \ip). 
Phase information is read out by applying a conditional phase 
shift J or I and interfering the paths on a beam splitter with 
a suitably chosen transmission (c) Experimental realiza- 
tion of the linear network. The input to the setup are either 
ultraviolet 390 nm (blue solid line) or fundamental 780 nm 
(dashed red line) pulses. Ellipses mark locations of effective 
beam splitters acting on orthogonal polarizations. For details, 
see Methods. 



quantum estimation theory [201 ] . the precision 5<p of the 
estimate is restricted by the quantum Cramer-Rao bound 
> 1/Vf, where F is the quantum Fisher information. 
In the case of a phase measurement, F is proportional to 
the variance of the photon number fi s sent through the 
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sample: 

F(|^)=4[^|^|^-((^|n s ^)) 2 ]. (1) 

Thus phase sensitivity is intimately related to the un- 
certainty of the photon number in the sensing arm. 
If the total number of photons is limited to N, this 
variance is maximized by a NOON state @, H| of the 
form (]N0) - \0N))/V2, describing a coherent, equally 
weighted superposition of all the photons traveling in 
either the sensing or the reference arm of the interfer- 
ometer. Here kets enclose a pair of figures that specify 
the numbers of photons in the two interferometer arms. 
An application of NOON states improves the precision to 
1/N, known as the Heisenberg limit [ll[, which defines 
the ultimate lower bound for phase estimation [2l|. This 
advantage however fades away if the phase shift is ac- 
companied by attenuation: even loss of only one photon 
destroys completely the superposition and the associated 
capability to acquire phase information [lH , [Tij . 

The simplest non-trivial example of a NOON state is a 
superposition ( 1 20) — |02))/\/2 which can be prepared us- 
ing the Hong-Ou-Mandel effect [22[ , by sending two pho- 
tons onto a balanced beam splitter. After transmission 
through the phase element the two superposition terms 
acquire a relative phase 2gg, leading to a doubled density 
of the interference fringes [23| which is behind the preci- 
sion enhancement. When a photon is lost in the sample, 
the quantum state of light inside the interferometer col- 
lapses to 1 10), which does not have any phase sensitivity. 
In order to avoid this problem, we take the probe state 
in the form of a three-component superposition 

|V>)=V^|20) + ^T|11)-V^|02), (2) 

where the non- negative weights xq, xi, and xi sum up to 
one. Let r\ be the intensity transmission of the sample. 
We will denote by p\ the probability of losing / photons 
and write the corresponding conditional state as ^/pil^i), 
where the ket \ipi) is normalized. When zero or one pho- 
ton is lost, the conditional states read: 

V^o#o) = V^|20) + v^HII) - V^|02) 
V^TlVi) - V2r?(l - V )X2\10) + V(l - r?)xi|01) (3) 

while only for the loss of both the photons the phase 
sensitivity is destroyed. Because the conditional states 
are distinguishable by the total photon number, the 
quantum Fisher information is given by a combination 
p F(\ip )) +p 1 F(\ip 1 }) [lj|. Maximization of this quan- 
tity yields optimal values of xo, x%, and X2, depicted in 
Fig. [2ja). The resulting precision is compared with the 
SIL and that achievable with original two-photon NOON 
states in Fig. [2](b), clearly showing preservation of the 
non-classical enhancement. 

In order to demonstrate the robustness of the opti- 
mal probe states we built a passive linear-optics net- 
work shown schematically in Fig. [TJb), fed with photon 
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FIG. 2: (a) Weights X2, x\, and xo of components in the op- 
timal two-photon probe state defined in Eq. ([5} as a function 
of the phase shift transmission r\, shown together with the 
probability p of preparing the optimal state using the scheme 
shown in Fig. \V[b) based on two beam splitters; (b) the pre- 
cision of phase estimation using two-photon optimal states 
(solid blue), NOON states (dashed red), and the standard in- 
terferometric limit (dotted gray) . The shaded area depicts the 
non-classical region, and the four arrows mark transmissions 
rj for which experimental data have been taken. 



pairs 1 11) produced via spontaneous parametric down- 
conversion 24]. Preparation of the optimal state is ac- 
complished by sending a pair of photons to a beam split- 
ter characterized by an intensity transmission which 
produces a state y / 2i} 1 (l - ^) (|20)-|02))+(2i? 1 -l)|ll). 
In particular, taking i?i = , yields the two-photon 
NOON state. In order to change the weights between 
the components 1 20) and |02), a second beam splitter 
with an intensity transmission $2 is inserted into one 
of the emerging beams. Provided that both the pho- 
tons are transmitted, the prepared state has the unnor- 
malized form ■ v /2i9 1 (l — i?i)|20> + V^(2^i - 1)|11) - 
^2\/ 2 ' t5 'i( 1 — ^i)|02). A suitable choice of i? x and d 2 
allows us to prepare optimal states with the success 
rate shown in Fig. [2a). Next, one of the produced 
beams undergoes a phase shift ip accompanied by a con- 
trolled loss introduced by the non-unit transmission rj 
of a beam splitter. Lost photons are monitored with a 
photon counter C. The attenuated state enters the detec- 
tion stage, which is designed to emulate a measurement 
saturating the Cramer-Rao bound at = 0, thus opti- 
mally extracting the encoded phase shift. The two beams 
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FIG. 3: Two-fold coincidence rates measured with photons 
prepared in (a) the optimal state and (b) the NOON state for 
the transmission r\ = 0.361 and 15 phases around (p = 0, 
taken with either ^ (coincidences AC, BC, CC) or -| (coinci- 
dences AA, BB, AB) additional phase shift before the beams 
interfere. The curves depict fitting of a model including ex- 
perimental imperfections. Vertical arrows mark phase values 
used to obtain histograms in Fig. [4] 



are interfered on a beam splitter with an intensity trans- 
mission followed by photon counters A and B. If 
the probe state does not include the |11) component, or 
only one photon makes it to the detection stage, we take 
i?d = 5, while if neither of these two conditions holds, the 
beam splitter transmission needs to be adjusted for the 
specific value of losses. Moreover, reaching the Cramer- 
Rao bound requires a conditional phase shift before the 
beams are brought to interfere: ? if no photon is lost 
and ~ if only one photon makes it through the interfer- 
ometer. This adjustment shifts the steepest slopes of the 
coincidence fringes to the vicinity of tp = 0. The adaptive 
nature of the optimal measurement requires in principle 
a joint quantum nondemolition photon number measure- 
ment [25j on both the beams after probing the phase 
shift, and then adjusting the beam splitter characteris- 
tics. In a passive linear optics implementation, this may 
be simulated by collecting separate data for two settings 
and postselecting events according to the number of pho- 
tons registered by the photon counter C. The limitations 
of the state preparation procedure, together with mode 
coupling and photon counting efficiencies mean that the 
measurements must be carried out in the coincidence ba- 
sis [26j]. The data used as the input for the estimation 
procedure are two-fold counting events for all the combi- 
nations of the photon counters. 

We implemented the linear optics network presented 
above using a pair of calcite displacers [I?) , as described 
in Methods. The two- fold coincidences are labeled with 
the paths taken by the photons at the output of the setup. 
Their dependence on the phase shift tp around the point 
of operation is illustrated in Fig. [3] for optimal and NOON 
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FIG. 4: Histograms of phase estimates for r\ — 0.361 obtained 
using photons prepared in (a) the optimal state and (b) the 
NOON state, probing two phase shifts differing by 0.12 rad. 
Each phase estimate was obtained from a series comprising 
approx. 2000 coincidences. 



states. The measurements were taken for four different 
values of loss corresponding to beam splitter transmis- 
sions of 7] = 0.2,0.361,0.4, and 0.547, employing two- 
photon optimal and NOON states, as well as using atten- 
uated laser pulses to establish the SIL. For each value of 
losses, we performed first scans over the entire 2ir range 
of the phase shift tp in order to calibrate the apparatus. 
The phase dependence of the coincidences counts was 
used to fit the parameters of the theoretical model [H| , 
which includes the imperfections of the state preparation 
and linear-optics manipulations. Next, precise measure- 
ments were carried out for 15 values of the phase around 
the most sensitive point ip = 0, with an increment of 
0.02 rad. In order to acquire information about the es- 
timation uncertainty, for each of the selected phases we 
took between 100 and 300 series of measurements, one 
series comprising on average 2000 coincidence events. 

A phase estimate was extracted from each series using 
the maximum-likelihood method 29], which is known to 
saturate the Cramer-Rao bound when applied to large 
samples. Histograms of phase estimates for two phases 
separated by 0.12 rad, obtained using experimentally 
generated optimal and NOON states, are shown in Fig. 0] 
along with fitted gaussian distributions. In order to gain 
a quantitative insight into the improvement, we depict in 
Fig.[5]the estimation uncertainty for five phases closest to 
(p = 0, when the sensitivity of the measurement scheme 
is highest. The uncertainties obtained from experimental 
data are rescaled by the square root of the number of reg- 
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FIG. 5: Uncertainties of phase estimates obtained using two- 
photon optimal (circles) and NOON (squares) states, as well 
as attenuated laser pulses in the SIL regime (diamonds), 
rescaled by the square root of the number of coincidences. 
For each transmission r/, data are shown for five phases 
ip — 0, ±0.2, ±0.4 rad. Horizontal lines represent the theo- 
retical Cramef-Rao bounds for given classes of input states 
taking into account imperfections of the interferometer. 



istered coincidences in order to obtain the figure of merit 
of the effective uncertainty per photon pair that can be 
compared with the theoretical Cramer-Rao bound. The 
advantage of using the optimal states is clearly seen for 
all the four values of losses used in the measurements. 
The results match well the theoretical model including 
experimental imperfections. 

Quantum-enhanced metrology requires new strategies 
to cope with decoherence, which are intrinsically dif- 
ferent from those targeting protection of optical qubits 
for quantum information processing applications [lg, [13] • 
Our experiment demonstrates the fundamental feasibility 
of such strategies. Further developments, such as exploit- 
ing multiphoton states [§, [l(| and going beyond the co- 
incidence basis (3(| , will allow a proper assessment of the 
capabilities of such approaches including all resources. 



Methods 

Photon pairs are generated in a 1 mm long bismuth bo- 
rate crystal X oriented for type-I process, pumped with a 
doubled output from a Ti:sapphire oscillator emitting a 
80 MHz train of a 780 nm central wavelength and 140 fs 
duration pulses. The photons, collected from the oppo- 
site ends of a cone with the half-opening angle of 4° are 
converted into orthogonal polarizations using a half- wave 



plate HWP, combined on a polarizing beam splitter PBS, 
filtered through a 5 nm interference filter IF, and coupled 
into a single-mode fibre. The contribution of higher order 
photon numbers is negligible. The relative delay between 
the photons is controlled by a motorized trombone placed 
in the path of one of the photons. The fibre is wound on 
a manual polarization controller PC to recover the |11) 
state at the output of the fibre. We found an adequate de- 
scription for the the two-photon st ate a t the output of the 

e/2(|20) + |02», where 



fibre of the form y/T~ 

e is typically of the order of 0.05%, and the phase 5 is in- 
troduced by the fibre. The second term can be attributed 
to a combination of polarization-dependent loss and bire- 
fringence in the single mode fibre delivering photon pairs 
to the interferometer. Its presence is seen in phase scans 
for NOON states and non-zero losses shown in Fig.^b) as 
a faint modulation of coincidences AC and BC, which in 
principle should remain constant. The setup also enables 
to launch into the fibre 780 nm fundamental-wavelength 
laser pulses to determine the SIL. 

The actual implementation of the linear optics net- 
work, shown in Fig.fljc), is based on a pair of calcite dis- 
placers and a Clan- Taylor polarizer, which ensures high 
visibility and stability of interference. The beam splitters 
are realized in a common-path version operating on two 
orthogonal polarizations, with the splitting ratios defined 
by the orientation of half-wave plates located before the 
birefringent medium. The second displacer is mounted on 
a rotation stage with the axis perpendicular to the plane 
of the setup. The orientation of the displacer defines 
the phase shift to be determined. The photons leaving 
the setup are coupled into 50:50 multimode fibre cou- 
plers, whose ends are plugged into single photon count- 
ing modules. Coupling into multimode fibres is adjusted 
to equalize effective efficiencies of the modules. In this 
arrangement, only 50% of events AA, BB, and CC are 
detected. In order to compensate for this inefficiency, in 
the numerical postprocessing half of events AB, AC, and 
BC has been randomly removed. 

The indistinguishability of the interfering two-photon 
pathways has been tested by setting $i = \ and sending 
the two emerging beams to separate detectors. The ob- 
served depth of the Hong-Ou-Mandel dip [22j exceeded 
98%. For the subsequent measurements, the relative de- 
lay of the two photons was set at the dip minimum. 
The classical visibility of the interferometer was above 
98%. In order to model accurately experimental data, we 
developed generalized expressions for coincidence count 
probabilities that include effects of residual photon dis- 
tinguishability, non-unit interference visibility, and fibre 
delivery [H|. 
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